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ABSTRACT 



o 
O 



Aims. We model the deepest observable layers of dark sunspot umbral atmospjwres in terms of an empirical model 
which equally describes observed near infrared continuum intensities and line p^^los. 

Methods. We use the umbral continuum intensity at 1.67/i and the three CI lines at 16888A, 17449A and 17456A to 
model the deep layers near the minimum of H~ absorption. An extrapolation of -.ircbra models to such deep layers must 
assure that (a) the calculated 1.67/i continuum does not fall below the range ol-yt^servations and (b) the resulting CI 
lines do not come out stronger than the lower observational limit. 

Results. Our calculations show that a T^Itr) stratification yields the best compromise between both criteria: a flatter 
gradient violates (a) a steeper one (b). We determine T^ff from the umbral and pBSjtospheric flux ratio down-scaling the 
monochromatic photospheric flux with the umbral contrast for each frequency v aS24)btain the monocharomatic umbral 
flux. Integration of both over u gives the ratio of total umbral and photospheric flujL which yields 3560 < T^f f < 3780 K. 

We assume for our model Te// = 3750 K and fit it to the theoretical model by ]\jgyer et al. (1974). Comparison of the 
gradient V = {dlogT) / {dlogPg) with the adiabatic one shows that umbral conveietion, if existing at all, can only occur 
at considerably deeper layers than in the photosphere. ^ 
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1. Introduction 

Recent observations of umbral continuum contrasts in the 
near infrared, carried out by Hall (1970) and by Ekmann 
and Maltby (1974), have led to a re-discussion of the 
temperature and pressure stratification in sunspot umbrae 
(Zwaan, 1974; Kjcldseth-Moe and Maltby, 1974). Their 
models disagree considerably with the empirical model MO 
by Stellmacher and Wiehr (1970, hereafter referred to as 
Paper I), which was shown to optimally fit the profiles of 
several magnetically insensitive lines observed in umbrae 
and had been confirmed by Kneer (1972) and by Koppen 
(1974). In Stellmacher and Wiehr (1971, hereafter referred 
to as Paper II) this fit was extended to the center-to- limb 
variation of a non-split line, which represents a very strong 
criterion for empirical models. 

MO had been deduced from line profiles and continuum 
contrasts at 5000A< A < 8000A, and can thus not rep- 
resent the deep layers ro.a > 1.62, for which it gives an 
arbitrary extrapolation. The present study demonstrates 
how MO can be extended towards deeper layers to fit ob- 
served IR data. The still existing scatter of observational 
data and of the photospheric models (to which any umbra 
model necessarily refers) allows an extention of MO to deep 
layers preserving the good line profile presentation in pa- 
pers 1 and 11 without assumption of additional parameters. 
We consider our empirical model as a 'boundary condition 
for three dimensional sunspot models' and for the 'quanti- 
tative interpretation of line profiles' (Zwaan 1974). 
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2. Observationa|1n;lata 

a\ 

Considering the 'S^chal decrease' of observed umbra-to- 
photosphere intensities in the visible region from 0.24 (at 
5452 A, Michard(353) to 0.04 (at 5790 A, Ekmann and 
Maltby 1974), recently observed infrared data may hardly 
be considered as ffidal results'. In particular, the determi- 
nation of an empTn^al umbra model from the low data in 
the visible range ^^ether with recent high infrared data 
will be uncorrectr.'^esides, differences between individual 
umbrae require tS^Jrse of data from one spot taken under 
identical conrfiiiorfasuch as the simultaneous broad band 
measurements by Ekmann and Maltby (1974). 



Table 1. Correction factors for absorption lines in the 
broad-band wavelength regions used by Ekmann & Maltby 
(1974). 



A-region 




0.576-0.582 


1.050 


0,664- 0.674 


1.092 


0.866-0.886 


1.00 


1.195-1.235 


im 


L630- 1.710 


1.044 


2-100-2300 


1.021 
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We correct their data for the influence of Fraunhofer 
lines integrating the ratio I^^brajjphot ^^^^ jj^jpg (■;l974) 
infrared atlas in the respective A- windows. For the visible 
region the same is done for the data by Wohl (1975). The 
correction factors, given in Table 1, might still be underes- 
timated due to the uncertain continuum level. 




Fig. 1. Observed relative umbral continuum intensities 
compared to calculations with our former model MO [dashed 
line) and with the present model M3 (solid, line, both rela- 
tive to Labs & Neckels 1968); dash points = M3 relative to 
Peteyraux (1952) and Badinov et al. (1965). Bars give raw 
and maximum corrected observations by various observers. 



We finally applied these corrections to the data from all 
umbrae at ^? < 20° given by Ekmann and Maltby (1974). 
The results are given in Fig. 1 together the umbral intensi- 
ties by Hall (1974), which appear to be considerably higher. 
Hall's data refer to continuum windows and are thus not 
affected by Fraunhofer lines; only parasitic light may inter- 
fere which docs not exceed 2% considering the Fell lines in 
that A-region to be of non-umbral origin (cf., paper I). The 
remaining differences to the Ekmann and Maltby (1974) 
data cannot be analyzed, since Hall (1974) did not observe 
in the visible A-range; light bridges in his spots may also 
explain his higher intensity values. 



3. Extrapolation of MO to the improved M3 model 

It has been shown by Zwaan (1974) that the strong varia- 
tion of relative umbra intensity between the maximum H~ 
absorption at 0.8 fi and the absorption minimum at 1.6 /x re- 
quires the steepest possible temperature gradient, i.e. radia- 
tive equilibrium whith = 3/4 • Tfjj{TR + q)-, where Tpj j 
and T/j denote the effective temperature and the Rosseland 
opacity, respectively and q is constant for deep layers. We 
therefore assume a T^{Tii) gradient for the extrapolation of 
our model MO to deep layers. 

To estimate the actual umbral T^f f value we multiplied 
the monochromatic photospheric flux Fp'^°* (Unsold, 1956, 
Fig. 20a) by the observed monochromatic intensity ratio 
jumbra I jphot ^ This yields the monochromatic umbral flux 
pumbra ^^der the assumption of a vanishing center-to-limb 



variation of the contrast at all i^, which is found in ob- 
servations and in all umbral models under discussion (cf., 
paper II). Planimetry over the frequency u then yields for 
umbra and photosphere their total (integrated) flux F^'^^"' 
and FP^°*. 

Assuming now radiative equilibrium, r.e., it is 

pumbra I pphot ^ T^jf{umhra)/T^ff{phot), and we obtain 

3560 K< Tg"™*-™ < 3780 K; the scatter reflects the range of 
observed contrasts. In order to keep the temperature gra- 
dient as high as possible we assume J'^^™'"'" = 3750 K; the 
resulting model M3 is listed in Table 3, it is shown in Fig. 2 
together with MO and some other models. The validity of 
this M3 has now to be checked by comparison with line 
profile and continuum observations. 
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Fig. 2. Umbral models: Zw=Zwaan (1974), 
KMMd=Kjeldseth-Moe & Maltby (1974) dark, MO 
our former, M3 the present model, 'fit' = adaption to 
Meyer et al. (1974); compared to photosphere models: 
Holw=Holweger (667), HSRA=Gingerich et al. (1971). 



3.1. The M3 representation of the continuum 

In order to avoid uncertainties from the deep layers of pho- 
tospheric niodcils we normalize our calculated umbra inten- 
sities exclusively to observed photospheric absolute inten- 
sities: Labs and Neckel (1968) list for 1.24/x < A < 3.6// 
the relative measurements by Pierce (1954), which they 
calibrated using the 'model distribution between 1.0 /U and 
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1.5 /i as a linking medium'. This affects that the two inten- 
sity values given by Pierce for 0.99 fi and for 1.05 /x are 6% 
too high. If we accordingly correct the relative intensities 
by Pierce (1952), we obtain a remarkable agreement with 
Peyturaux (1952). 

A similar uncertainty concerns the photospheric ab- 
solute intensities at 3.6/i < A < 4.8/i which, according 
to Koutchmy and Peyturaux (1970) show a large scatter 
among the various observations. In particular, Badinov et 
al. (1965) publish absolute intensities, which are consider- 
ably lower than those listed by Labs and Neckel (1968). 

Consequently wo normalize for those two wavelength 
regions the intensities from our umbra model alternately to 
both, the data listed by Labs and Neckel (1965) and the 
lower values given by Peyturaux (1952) and by Badinov 
et al. (1965), respectively. The resulting relative umbral 
intensities for MO and for M3, given in Fig. 1 together with 
the observations discussed in Sec. 2, show that model M3 
fully agrees with the range of continuum observations when 
considering those uncertainties in the photospheric data. 

3.2. The M3 representation of line profiles 

Since model MS differs from MO only in the deeper layers 
To. 5 > 1.6 (see Fig. 2 and Table 3) one would not expect 
significant influence on the line profiles mentioned in Sec. 1. 
However, calculations show that the contribution of these 
deeper layers to the line wings of NaD2 and of Fe 5434 is not 
sufficiently small to be neglected. Moreover, the line profile 
calculations show that M3 represents an 'upper limit' for 
an extrapolation of MO to deep layers. 

For our line profile calculations we adapted the damping 
constants 'y/jt'^^"'^ to the observed photospheric line wings 
taking the known (g ■ /) values by Lambert and Warner 
(1968) for NaD2 and by Wolnik et al. (1970) for Fe5434 to- 
gether with the relative abundances log{eNa/eH) = —5.75 
and log{eFe/£H) = —4.7, thus improving our fit of log{g f e) 
to W^^°* in papers I and II. With 7jvaD2 = 1.4 7*^^^°'' and 
7Fe5434 = 1.87*'"^°'' wc achieve an agreement within 1% 
with the observed photospheric line wings at AAFe5434 > 
0.1 A and at A.XNaD2 > 0.3 A, when using Holweger's 
(1967) photospheric model and assuming LTE. 

The validity of LTE is reasonably be also assumed for 
the umbral model. The line profiles obtained with M3 are 
shown in Figs. 7 and 8 together with our observations dis- 
cussed in Papers I and II. The bars give the raw data as 
upper and the maximum uncorrected ones as lower end. 
The maximum correction corresponds to actually vanish- 
ing Fe II 6149 and Fe II 7224 lines in the dark umbral core. 
Figs. 7 and 8 show that the models by Zwaan (1974) and 
by Kjeldseth-Moe and Maltby (1974) yield profiles outside 
that range. 

As a further criterion for the temperature gradient 
in deeper umbral layers we use the equivalent widths of 
C 1 16888A, C 1 17456A and C 1 17449A which we determine 
from Hall's (1974) infrared atlas to 7.5, 13.0, and 20.0 mA. 
These values represent upper limits since any parasitic light 
would have strengthened them. The qiiantity (g / e) was de- 
termined for the C I lines by a fit to the photospheric equiv- 
alent widths from Hall's atlas: 14^^''°* (C 16888) = 94 mA, 

wf°\C17456) = 134mA and W|''"'*(C 17449) = 168mA 
(see Table 2). 



For the calculation of the C I lines one has to consider 
possible influences of the formation of CO on the num- 
ber density of atomic carbon. That of CO depends slightly 
on the formation of OH and H2; the influence of other 
molecules is negligible. For M3 we flnd nco ~ "-c at the 
To. 5 ~ 5.0 layer, which also gives the maximum contri- 
bution to the C I lines and their adjacent continuum. The 
eventual influence of CO-formation on the C-lines amounts 
to 25%. The resulting VT^™''™ (Table 2) indicate that model 
M3 represents the 'steepest possible extrapolation' of MO 
and gives the best fit to the observations. 




Fig. 3. Line blanketing and deviations from rad.equ. for 
the present model M3 and the photospheric by Holweger 
(1967). 




Fig. 4. Depth dependence of Rosseland opacity ratio re- 
quired to convert the blanketed into the unblanketed mod- 
els. 
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Fig. 5. Fit of model M3 to the theoretical umbra stratifi- 
cation by Meyer et al (1974) without affecting the 1.67yU 
intensities. 



4. Conclusion 

Having shown that M3 well represents observations of line 
profiles and continuum, it is reasonable to discuss that 
model in more detail, keeping in mind that the validity 
of M3 is limited by the accuracy of the observations the 
model is based on. 

In order to obtain some idea about the umbral line 
blanketing we compare M3 with the corresponding unblan- 
keted gray r.e. stratification with Tg// = 3750° K (cf., 
end of Sec. 3). A presentation in the y4j'^contjn-j pja^Q 
(Fig. 4) shows that M3 has a typical line blanketing with 
'back-warming' at 1.65 > tr > 10~^ (corresponding to 
2-3 > To. 5 > 3 • 10""^) and 'lowering of the boimdary tem- 
perature' at tr < 10^'', as is predicted by the picket- fence 
model (Chandrasekhar, 1935). 

A more quantitative investigation of the line blanket- 
ing can be achieved by transforming the blanketed into the 
unblankcted model. The resulting quantity k/j/k^"* (Fig. 
5) shows for M3 a much stronger depth dependence than 
for the photosphere. Such a behavior has already been dis- 
cussed by Mattig and Schroter (1974). 

An interesting question is the limit of radiative equilib- 
rium (r.e.) below which convection is possible. For a rough 
estimate we connect M3 to the theoretical model by Meyer 
et al. (1974) in the 6 = 5040/T versus logPg plane (Fig. 6). 
This is done keeping the 1.67/i intensities unaffected within 
10^"^ and requires a validity of r.e. up to T0.5 < 8. The ex- 
isting uncertainties in the observed photospheric absolute 
and umbral relative intensities do not allow us to exclude 
an onset of deviations from r.e. at slightly higher layers, 
which might then give rise to observable umbral convection 
near 1.67/Lt. 

Another criterion for convection is the maximum in the 
depth dependence of V = {dlogT)/{dlogPg) (Fig. 7). For 
the umbral M3 model the V-maximum occurs as deep as 

^umbra 
~0.5 

phot ^ 
^0.5 



a 10.0 but for the photosphere much higher at 
"^0.5 ~ 2.0. On the other hand, the criterion for con- 
vective instability, V = Wad = 0.4 is satisfied for M3 at 
^umbra ^ 2.4 and for the photosphere at rff = 0.44 (see 
Fig. 6). 

Both considerations indicate that convection is possible 

for umbrae only at considerably larger depths than for the 
photosphere, in full agreement with the absence of macro- 




Fig. 6. Depth dependence V indicating the much deeper 
onset of convection in the umbra as in the photosphere. 

turbulence in line profile representations with the new em- 
pirical working model M3 for sunspot umbrae. 
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. 7. Fe5434 at cos'd = 1.0 (upper) and cos^? = 0.5 {lower panel); calculated {solid lines); observations {bars). 
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Fig. 8. Calculated and observed profiles of the violet wing of NaD2; x = raw data by J. Harvey (priv. com.) 
Table 2. Umbral CI lines from observations {upper) and from various emirical models {lower part). 



C 16888.3 



C 1745&.0 



C 17448.8 



W^"' (Hall atlas) 
tyumbra (Hall atlas) 



94. 
- 3.585 
7.5 



134. 
- 3355 
13. 



168. mA 
- 3.170 
20. mA 



Calculations with umbra models: 

M3-fiiial 7.2 

M3-r* 7.9 

MO a3 

KMMd 0.6 

KMMb 3.0 

Zw 13.3 



11.2 
12.4 
0.4 
1.0 

20.1 



15.9 
17.5 

a6 

1.4 

6.6 

27.8 
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Table 3. Model M3 with deep layers as stratification (upper) and fit to the Meyer et al. model for r > 1.0 (lower 
part); also given are the abundances used. 
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